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The invention pertains generally to data transmission and reception, and, more particularly, to the reception 
of high speed data that has been used to modulate many carriers. This method of modulation is know by many 
names-Orthogonally Multiplexed Quadrature Amplitude Modulation (OMQAM), Dynamically Assigned Multiple 
QAM (DAMQAM), Orthogonal Frequency Division Multiplexing (OFDM); herein it is referred to as multicarrier 
5 modulation. 

In general, oscillators (frequency sources) in the transmitting and receiving modems of a data communi- 
cation link cannot by synchronized, and the receiving modem must use special circuitry or signal processing 
algorithms to of the received signal; this task is often referred to as clock and carrier recovery. 

The task is complicated by the introduction in some transmission media, most notably the General Switched 

10 Telephone Network (GSTN) ( of frequency offset and phase jitter as shown in Fig. 1. Frequency offset, desig- 
nated herein as f, is the difference between the frequency(ies) of the received carriers) and the reference car- 
riers) generated in the local receiver it may be caused by (a) a frequency difference between the carrier(s) 
used in the transmitter and the reference carriers) in the receiver and/or (b) a mismatch between the fre- 
quencies of the modulating and demodulating carriers in Frequency Division Multiplexing (FDM) equipment in 

15 the network; on the GSTN the resultant combined offset may be as much as 5 Hz. Phase jitter is seen as phase 
modulation of the received signal, and often has a few discrete, identifiable frequency components; the power 
supply frequency and telephone ringing frequency (60 Hz and 20 Hz, respectively, in the U.S.A.) are common 
components. The carrier recover circuitry or algorithms must track, or follow frequency offset and phase jitter 
in order to ameliorate their effect on the received signal. 

20 Circuitry and algorithms for recovering clock and carrier (including tracking frequency offset and phase jit- 

ter) from data-modulated single-carrier signals have been well documented in the literature (see J.A.C. Bin- 
gham, The Theory and Practice of Modem Design, John Wiley, 1988). Nearly all of the prior art uses 
phase-locked loops (PLLs) of one form or another, and designers have long recognized the harmful effedts of 
delay inside those loops. For carrier recovery it was agreed that tracking of phase jitter is made much more 

25 difficult, maybe even impossible, if the delay of an adaptive equalizer (typically about 10 ms) is included in the 
loop; therefore special algorithms (see, for example, D.D. Falconer, "Jointly Adaptive Equalization and Carrier 
Recovery in Two-Dimensional Communications Systems", BellSvst Tech J., vol. 55, pp. 317-314, March 1976) 
were developed to remove this delay from the loop. 

In multicamer modulation, data are grouped into blocks of bits; in the systems described by Baran (U.S. 

30 Patent 4,438,511) and Hughes-Hartogs (U.S. Patent 4,679,227) the blocks may comprise more than one 
thousand bits. Each carrier is modulated by just a few of those bits, and the modulation is held constant for the 
duration of one block; this duration, or symbol period, may therefore be several hundred times the symbol period 
of a single-carrier modem. Furthermore, the symbol period may be much greater than the periods of the com- 
ponents of the phase jitter that are to be tracked. 

35 Signal processing in the multicamer receiver must be performed at the symbol rate, and information about 

the received signal-data contained therein and imperfections (frequency offset, phase jitter, etc.) thereof-is 
available only after each block has been performed at the symbol rate, and information about the received sig- 
nal-data contained therein and imperfections (frequency offset, phase jitter, etc.) thereof-is available only after 
each block has been processed. This delay of one symbol period, which may be as large as 130 ms (more 

40 than ten times the delay through most adaptive equalizers), therefore appears inside the carrier recovery loop, 
and makes conventional jitter tracking impossible. 

The long symbol period associated with multicarrier modulation makes the tracking of phase jitter imposed 
on multi-carrier signals a much more difficult problem than for single-carrier signals. One proposal (B. Hirosaki, 
et al, "A 19.2 kbit/s Voiceband Data Modem Based on Orthogonally Multiplexed QAM Techniques", IEEE Intl. 

45 Conf. Commun. Rec, pp. 661-665, Aug. 1985) was to input an unmodulated pilot tone to a set (one for each 
carrier) of adaptive jitter predictors, and feed the output signals forward to cancel the jitter on each modulated 
carrier. This approach has several disadvantages: the amount of information about the jitter available from one 
pilot placed at the edge of the available frequency band is very small, the tapped-delay-line form of the predic- 
tors is poorly suited to filtering single-tones, and the approach requires a large amount of computation-par- 

50 ticularly if a large number of carriers is used. 

Another problem is that the jitter frequencies are not Subspace Approach to Multiple Emitter Location and 
Spectral Estimation", Ph.D. Thesis, Stanford University, 1981 , and the ESPRIT algorithm described by R. Roy, 
A. Paulraj and T. Kailath in "ESPRIT-A Subspace Rotation Approach to Estimation of Parameters of Cisoids 
in Noise", IEEE Trans. ASSP . vol. ASSP-34, p. 1 340-1 342, Oct 1 986) but they all require a very large amount 

55 of computation, and the particular problem of determining the frequencies of phase jitter has not been addres- 
sed. 

Many of the early objections to multicarrier modulation were based on the assumption that correction for 
phase jitter was impossible because the problems of long symbol time and unknown jitter frequencies had not 
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been, and probably could not be, satisfactorily solved. 

In accordance with the preferred embodiment of the present invention, the problem of the symbol delay 
occurring inside the earner recovery loop is solved by combining the concepts of the block processing needed 
for the efficient detection and decoding of multi-carrier signals with those of the serial processing needed to 
5 generate the local correcting signals that have the required compensating offset and jitter. In the present inven- 
tion this is done by using feed-back (also known as closed-loop) algorithms or circuits in a combination of serial 
and block modes. 

These feedback algorithms are designed, however, primarily for continuous operation while receiving da- 
ta-modulated signals; their adaptation is fairly slow, and they presuppose a knowledge of the frequencies of 
w the phase jitter components that are to be tracked. 

Therefore, upon connection of a channel forthe transmission of data, the data signal is preceded by a train- 
ing signal comprising a few pilot tones. In the receiver these tones are separated by filtering, their phases are 
calculated as functions of time, and those functions are analyzed to form initial estimates of (a) dock offset, 
(b) carrier offset, and (c) the amplitudes, frequencies and phases of all significant components of phase jitter. 
is These estimates are then used to initialize the circuits or algorithms which generate a phase-correcting 

signal. This signal is used during the reception of data-modulated signals as follows: 

1. In the serial mode, samples of the received signal are multiplied by the phase-correcting signal, which 
comprises components at the jitter and offset frequencies, to produce samples of a corrected signal. 

2. A block of these corrected samples is then demodulated and decoded (as described, for example, in 
20 U.S. Patent 4,679,227) to produce a block of estimated received data. These data, togetherwith the demod- 
ulated (but not decoded samples), are used both to update the estimates of the amplitude and phase 
characteristics of the channel, and as inputs to the algorithm or circuitfor estimating carrier and clock offsets 
described in paragraph 5 below. 

3. The block of received data is then remodulated to produce a sequence of samples of a reference signal, 
25 and in a serial mode the phases of these samples are compared to the phases of the corrected samples 

to generate a sequence of samples of phase error. 

4. These samples of phase error are then analyzed to calculate required updates for the parameters of the 
algorithm or circuit that is used to generate the phase-correcting signal. 

5. The differences between the demodulated and the decoded data are analyzed and filtered to form esti- 
30 mates of the remanent carrier and dock frequency offsets. The estimate of carrier offset is used to update 

a parameter of the circuit or algorithm that generates the phase-correcting signal and similarly, the estimate 
of dock offset is used to control either the original sampling time of the received signal or an interpolator 
which operates upon the samples. 

6. The phase-correcting signal and sampling time thus updated are then used for reception of the next block 
35 of data signals, and steps 1 through 5, described above, are repeated for each subsequent block of data 

signals. 

For a better understanding of the present invention and to show how the same may be carried into effect, 
reference w9l now be made by way of example, to the accompanying drawings in which: 

Fig. 1 is a schematic block diagram representation of a data transmission link subjected to clock and carrier 
40 frequency offsets and carrier phase jitter. 

Fig. 2 is an overall block diagram of an embodiment of the present invention showing the feed-forward train- 
ing mode and the feed-back data mode. 

Fig. 3 shows a typical phase/time plot for the two pilot tones for the case of just one jitter component 

Fig. 4 shows a look-up table implementation of a generator of one sinusoidal component of the phase-cor- 
45 recting signal, with means for updating the amplitude and phase. 

Fig. 5 shows a typical phase/frequency plot for many carriers. 

The transmitter and receiver of a data transmission link, such as is shown in Fig. 1, each indude an oscil- 
lator (usually crystal-controlled) to control the frequency at which the transmit and receive signals are sampled, 
the data rates of those signals, and the frequencies of the carriers that are used for modulation. 
so Circuitry or signal processing algorithms are used in the receiver to lock the sampling rate in the receiver 

to that in the transmitter. Clock offset, e, is defined by the ratio of the natural (i.e., unlocked) sampling rate in 
the receiver, f w , to that used in the transmitter, f 8 : 

f„ = (1 + e)f. 0) 

Circuitry or algorithms are also used to lock the frequencies of the carriers used for demodulation to those 
55 inherent in the received signal. The common frequency difference (offset) between the received carriers and 
the natural, unlocked frequencies of the receiver, which may be caused by a combination of the differences 
between the frequencies of (a) the transmitter and receiver osdllators. and (b) the oscillators in the modulator 
and demodulator of the FDM equipment (shown as f«„ and fcJ, is defined as Af. 

3 
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The oscillators in the FDM equipment may be subjected to phase jitter; the sources are typically additive 
and their sum affects the phase of the modulating and/or demodulating carrier as shown in Fig. 1 . To simplify 
the explanation of the invention, it will be assumed that there is significant phase jitter at only two frequencies, 
ft and f 2 , and that the jitter components have amplitudes A t and A 2 , and phases 0 n and G 2 . Each component 
5 of the received signal at a frequency f, could then be represented: 
without jitten 

sin[a>t + 0] 

with jitter 

to sin[©t + 0 + A^sin^t + e^) + A^in^t + €>&)] 

where, for all frequencies (of signal and jitter), 

co = 2rrf 

As shown in Fig. 2, a real received signal, x(t), which may have been corrupted in this way, is first sampled 
by switch 10 at a rate f ar (= 1/T r ) and with a timing offset, t, and converted to a complex signal, Xp(nT r +x) + 
is jx q (nT r +r), by a Hilbert transformer 20 (see, for example, LR. Rabiner and R.W. Schafer, "On the Behavior of 
Minmax FIR Digital Hilbert Transformers", Bell Syst. Tech. J. t vol. 53, pp. 363-390, Feb. 1974). Upon first con- 
nection of a data channel the training algorithm is used; this is represented symbolically in Fig. 2 by switch 30 
in the Training - position. 



20 Training Mode 

The training signal, from the transmitter of Fig. 1, comprises two pilot tones at frequencies U and f 2 (1000 
and 2000 Hz in the first embodiment). It is sampled by switch 10, transformed by the Hilbert transformef 20, 
and then processed as follows: 
25 1 . The real and imaginary parts (also shown as the in-phase and quadrature components) of these tones 

are separated by band-pass filters 40-70 to generate the sampled signals y pl (n), y q i(n), y^n), and y^n). 

Ideally the response of each of these filters should be maximum at one of the pilot tone frequencies and 

zero at the other. Typical transfer functions-in terms of the z-transfocm variable, z, are 



30 



35 



F t (2> = 



and 



X p (2) 



2T - 2COS9 2 2 + 1 
2* - 2X0030,2 + r c 



(2a) 



40 



F 2 (2) 



XP(2) 



X q (2) 



z z - 2rcos9 1 z + l 
z z - 2rcose 2 z + 



(2b) 



where 0( = 2*f|/f s for i = 1 , 2. 

45 2. The output signals of the four filters are applied to arctangent operators 75 and 80, which calculate the 

phases of each pair of samples, taken at times t = nT r +r, the two separated tones according to 

<fc(n) = arctan^nj/yp^n)); (3) 
The phases of a number (preferably an integer power of 2, and equal to 1 024 in the preferred embo- 
diment) of samples are calculated sequentially and stored; they are then processed as a block in processor 
so 85 as described in steps 3-5. 

3. Each phase, <|>i(n) and (^(n), is fitted to a linear function of time (measured in increments of T r by the 
variable n) using the well known Least Mean Squared Error (LMSE) criterion (see, for example, B. Widrow 
and S.D. Steams, Adaptive Signal Processing, Prentice Hall, 1985); that is, 

<Mn)«a 1n + Mori = 1,2, (4) 
55 where a lf p t , a 2 , and p2 are the LMS estimates of the slope and zero-time intercept of each phase. An 

example of phases and <fc and the best linear fit to them is shown in Fig. 3 where, for sake on clarity, 
only one jitter component is shown. 

4. The carrier and dock frequency offsets, A f and e, are then estimated from the slopes, and a 2 : 
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a,f z - a 2 f t 

At = (5a) 



and 



2ir (f 2 - f,) 

10 l + € (5b) 

f, (a 2 - a,) 

5. A spectral analysis (a Fast Fourier Transform (FFT) was used in the first embodiment) is then performed 
is on the residual, non-linear parts of the phase functions, 

<fri'( n ) = <M") - ctin - Pi for i = 1 . 2 
to estimate the amplitude, A^ kt frequency, f Jk , and phase, ® lkl of the components of phase jitter (k =1, 2 for 
the illustrated case of two components). 

6. After a few blocks (the first embodiment uses just one) of the two pilot tones have been transmitted, 
20 received, and processed the estimated values of Aj k , f^, and © J>c are used to calculate the coefficients and 

initial states of the sinewave.generators 100 and 105. This completes the training mode. 
Switch 30 is then placed in the data position, and processing continues with the reception of data modulated 
onto multiple carriers. 

25 Data Mode 

Carrier Recovery: 

1 . Hilbert-transformed samples of the data-modulated multicanier signal, x(t) t which is assumed to be 
30 impaired by carrier frequency offset and phase jitter, are sequentially applied to multiplier 110, and there 

multiplied by samples of the correcting signal, cos$ c + jsin^c which are generated by the look-up table 260. 

The output signals of multiplier 110 (the corrected samples) are assembled into a block, which is 
stored in delay 270 and also processed as follows: 

2. The block of output signals from multiplier 1 10 is input to demodulator 120, which performs an FFT on 
35 the block of samples of the corrected signal, as described by E.O. Brigham in The FFT and Its Applications , 

Prentice Hall, 1974 (a general theory of FFTs), and by D. Hughes-Hartogs in U.S. Patent 4,679,227 (the 
application of FFTs to multicanier modulation). The output from 120 is applied to decoder 130, which esti- 
mates the transmitted data, as also described in U.S. Patent No. 4,679,227. 

3. The output signals from demodulator 120 and decoder 130 are applied to sub tractor 125; the output sig- 
40 nals from 125, which are the amplitude and phase differences between the demodulated and the decoded 

signals, are applied to processor 140, where a subroutine estimates the characteristics of the channel as 
described in U.S. Patent 4,679,227. 

4. The estimates from the sub-routine of processor 1 40. which are updated after the reception of each block 
of data, are then applied, together with the estimates of the transmitted data from decoder 130, to a remod- 

45 ulatbr 150 which performs an inverse FFT and generates a reference signal B. 

5. The phase differences from subtracter 125 are also applied to processor 1 60, where a linear fit subroutine 
calculates an LMSE fit of the phases to a linear function of frequency, af + b, as illustrated in Fig. 5. The 
constant term, b, is attributed to an error in the phase of the synthesized correcting signal, and a signal 
representative of b is applied to low-pass filter 1 70. 

so 6. Phase comparator 190 then calculates, sequentially, samples of a phase-error signal which is the dif- 

ference in phase between the corrected signal. A, which has been stored in delay 270, and the reference 
signal, B, from remodulator 150. If the two signals are represented as (x p (n) + jx^n)) and (y p (n) + jy<,(n)), 
respectively, then this phase difference is given exactly by 

55 
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A q (n)y p (n) - x p (n)y q (n) \ 

0_(n) = arctian / ) (8) 

( x p (n)y p (n) + x^n) y q (n) J 

However, because the loop was initializing during the training phase, the differences between the corrected 
and reference signals are small, and the arctangent function of equation (8) can be approximated by its 
argument 

Furthermore it was found that the division operation in the argument of the arctangent function can be dis- 
pensed with, and a simplified indicator of the phase-error, 

+e(n) « Xq(n)y p (n) - Xp(n)y q (n) (9) 
can be used; this is applied to multipliers 200 and 21 0. 

7. In multipliers 200 and 210 the samples of the phase-error signal, <|> ot are multiplied by samples of cosine 
and sine waves (i.e., complex samples) at the identified jitter frequencies, which are generated by look-up 
tables 90 and 95. The complex products of these multiplications are then applied to low-pass filters 220 
and 230, respectively. Filters 220 and 230 typically have transferfunctions of the form 

F(2) = (10) 

(Z - 1) (2 -Z 2 ) 

and the gain, zero, and pole (K, z 1( and Z2) are chosen according to conventional PLL theory, as described 
by F.M. Gardner in Phaselock Techniques, John Wiley, 1979, so as to optimize the transient response of 
the loop. 

The operations of multiplying and low-pass filtering together constitute a correlation of the phase error sig- 
nal with each of the complex jitter components. 

8. Steps 6 and 7 are performed for the duration of one block, and then the output signals of filters 220 and 
230 are sampled by switches 240 and 250, respectively, to produce an estimate of each component of the 
phase error of the correcting signal. These error components, which are calculated during one block will 
be used, however, to update the parameters of the correcting signal that will be used for the subsequent 
block. Therefore, because the block period is not necessarily an interger multiple of the period of either of 
the jitter components, each component must be shifted, in rotators 245 and 255, by an angle a> Jk T (=2Tcf jk T), 
which is the angle by which each identified jitter component will rotate in one block period. The sampled 
and rotated error components are then used to update the parameters of sinewave generators 100 and 
1 05, which then continue undisturbed during the next block. 

Sinewave generators 1 00 and 1 05 may be implemented in several different ways, and the preferred method 
of updating them depends on the implementation. In the preferred embodiment, shown in Fig. 4 the real and 
imaginary parts of the error components, e p and £q are converted to amplitude and phase errors, s A and s+, and 
these errors are then used to increment the amplitude, Aj k and phase, of each component of the correcting 
signal according to 

V = + 8e A k 

and 

where 5 is a small step-size multiplier. The updated phase, fa is then used as an address for a cosine and 
sinewave look-up table, and the two output signals are multiplied by the updated amplitude. 

9. The real output signals of sinewave generators 100 and 105 and of filter 170 are applied to a summer 
175 to form a composite correcting phase angle, 



2 

0 e =» 2r AfnT_ + 2 A Jk cos(2rf jk nT r + e jk ) 

which is the estimate of the phase shift required to compensate for offset and jitter. This <fr c is then applied to 
the sine and cosine look-up tables 260 in order to generate samples of a complex correcting signal. 
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Clock Recovery: 

The slope term, a, calculated by the linear fit algorithm 160, (described in Data Mode paragraph 5), is attri- 
buted to an error in the sampling phase, x, and is applied to low-pass filter 180. A typical transfer function of 
this filter, in terms of the z variable, is 



and the zero, z 1t and pole, z 2 , can be chosen according to standard PLL theory so as to achieve a compromise 
between transient response and noise bandwidth. 

The output signal of this filter may be used in one of two ways (as indicated by the dotted lines in Fig. 2): 

(a) it may be fed back to a Voltage Controlled Oscillator (VCO) 300, which is used to control the first sam- 
pling switch 10, or 

(b) the switch 1 0 may be controlled by an unadjusted (open-loop) oscillator, and the output of the filter inte- 
grated in integrator 290 and then used to control an interpolator 310 as described, for example, in R.W. 
Schaferand LR. Rabiner. "A Digital Signal Processing Approach to Interpolation", Proc., IEEE , June 1973. 
It will be apparent to those skilled in the art that the invention disclosed herein may be embodied in other 

specific forms without departing from the spirit and essential characteristics thereof. Accordingly, the disclosure 
of the present invention is intended to be illustrative, but not limiting of the scope of the invention that is set out 
in the following claims. 



Claims 

1. A method for estimating, and compensating for phase jitter imposed on a received multicarrier signal by 
imperfect components in a transmission medium, said method comprising the steps of: 

(a) multiplying said multicarrier signal received during a time interval by a locally generated correcting 
signal to produce a corrected received multicarrier signal, in which said locally generated correcting sig- 
nal is an estimate of the phase shift required to correct for the effects of said phase jitter; 

(b) generating a reference signal by demodulating and then decoding said corrected multicarrier signal 
to produce an estimate of the transmitted data, and then remoduiating said data to produce a reference 
signal; 

(c) calculating a phase-error signal which is the difference between the phase of said correcting mul- 
ticarrier signal and that of said reference signal; 

(d) filtering the calculated phase-error signal of step (b) to generate an estimate of the changes required 
to the parameters of said locally generated correcting signal, said parameters being the amplitudes, 
frequencies, and phases of a plurality of components of said correcting signal, or any equivalent set 
which contains the same information; 

(e) updating the parameters of step (d) by said changes to generate new parameters of said correcting 
signal; and 

(f) repeating steps (a) to (e) in each subsequent time interval. 

2. The method of claim 1 further including the steps of: 

(g) sampling said received multicarrier signal in synchronism with a locally generated clock signal; 

(h) assembling the samples of step (g) into a block; 

(i) performing the updating of step (e) after the reception of each block; and 

(j) multiplying the samples in a subsequent block of said received multicarrier signal by said correcting 
signal, which is defined by said new parameters. 

3. The method of claim 2 wherein said reference signal is calculated by the steps of. 

(k) demodulating and decoding said corrected received multicarrier signal to estimate what data was 
modulated onto each of the many carriers; 

(I) filtering the results of said demodulating and decoding in order to estimate the amplitude and cor- 
rected phase responses of said transmission medium at each of the many carrier frequencies; 
(m) remoduiating said carriers with said estimated data, using said amplitude and phase responses, to 
generate, as a reference, an ideal received multicarrier signal which has no phase jitter; and 
(n) comparing said corrected received multicarriersignal with said remodulated, reference signal to gen- 
erate said phase-error signal. 
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The method of claim 3 wherein said correcting signal comprises components at several frequencies, called 
jitter frequencies, at which it has been determined that there is phase jitter, further including the steps of: 
(o) correlating said phase-error signal with complex sinusoidal signals at said jitter frequencies; and 
(p) using the results of said correlations to define the changes that should be made to the parameters 
of several sinewave generators, which generate components, at the identified jitter frequencies, of said 
correcting signal. 

The method of claim 4 wherein said parameters are the addresses of look-up tables and the multipliers 
that are applied to the output signals of said tables in order to generate the components of said correcting 
signal. 

A method for estimating, and compensating for frequency offset imposed on a received muiticam'er signal 
by imperfect components in a transmission medium, said method comprising the steps of: 

(a) multiplying the multicarrier signal received during a time interval by a locally generated correcting 
signal to produce a corrected received multicarrier signal, in which said locally generated correcting sig- 
nal is an estimate of the phase shift required to correct for the effects of said frequency offset; 

(b) demodulating and then decoding said corrected received multicarrier signal in order to estimate what 
data was modulated onto each of the carriers; 

(c) filtering the results of step (b) in order to estimate the corrected phase response of said transmission 
medium at each of the carrier frequencies; 

(d) fitting said corrected phase response of the medium to a linear function of frequency; and 

(e) interpreting the zero-frequency intercept of said linear function as being caused by a difference be- 
tween the phase of said correcting signal and the phase that is common to all carriers of said multicahrier 
received signal jow-pass filtering and integrating said zero-frequency intercept, and using the result of 
said filtering to update the low-frequency component of said correcting signal. 

A method for correcting, in a modem that receives, via a transmission medium, data modulated onto many 
carriers, for a difference between the frequencies of a locally generated sampling clock and the clock in a 
remote, transmitting modem, said method comprising the steps oft 

(a) sampling said received multicarrier signal in synchronism with a locally generated clock signal; 

(b) interpolating between said samples according to an interpolation parameter in order to convert sam- 
ples taken at the frequency of said locally generated dock into interpolated samples at the frequency 
of the clock in said remote transmitting modem; 

(c) demodulating and then decoding said interpolated samples in order to estimate what data was mod- 
ulated onto each of the many carriers; 

(d) filtering the results of said demodulating and decoding in order to estimate the phase response of 
said transmission medium at each of the many carrier frequencies; 

(e) fitting said phase response of said transmission medium to a linear function of frequency; 

(f) interpreting the slope of said linear function of frequency as being caused by a difference between 
the sampling frequencies of transmitter and receiver; and 

(g) low-pass filtering said slope, and using the result as said interpolation parameter. 

The method of claim 7 in which said filtered slope is used as the control for a phase-locked loop which 
adjusts the frequency of said locally generated clock to match that of the remote transmitter clock. 

A method for calculating, in a modem that receives a multicarrier data signal which may be impaired by 
phase jitter, frequency offset and a difference between the frequencies of a locally generated sampling 
clock and the clock in a remote transmitting modem, the initial parameters of a locally generated correcting 
signal which corrects for said impairments, said method comprising the steps of: 

(a) calculating the initial values of the parameters of said correcting signal during the reception of a train- 
ing signal comprising a plurality of pilot tones; 

(b) separating said pilot tones by bandpass filters, and calculating their phases at successive sampling 
times; 

(c) analyzing said phases to identify linear and sinusoidal functions of time; 

(d) estimating the frequencies, amplitudes and phases of said sinusoidal functions of time; 

(e) using the estimates of step (d) to initialize the parameters of sinewave generators; 

(f) comparing the slopes of said linear functions of step (c), and generating an estimate of both said 
frequency offset and said clock frequency difference; 
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(g) using said estimate of said frequency offset to initialize the zero-frequency component of said cor- 
recting signal; 

(h) using said estimate of the clock frequency difference to initialize the state of a low-pass filter and 

(i) using the output signal of said low-pass filter to control an interpolator which interpolates between 
samples taken in synchronism with a locally generated clock signal to generate samples at the fre- 
quency of the dock in a remote transmitting modem. 

10. The method of claim 9 wherein the output signal of said low-pass filter is used to control a phase-locked 
loop which adjusts the frequency and phase of said locally generated clock to match that of the dock in a 
remote transmitting modem. 

11. The method of daim 1 wherein the samples of the receive signal are processed in blocks using Fast Fourier 
Transform (FFT) and inverse FFT (IFFT) algorithms. 

12. The method of daim 6 wherein the samples of the receive signal are processed in blocks using FFT and 
IFFT algorithms. 

13. The method of daim 7 wherein the samples of the receive signal are processed in blocks using FFT and 
IFFT algorithms. 

14. The method of daim 9 wherein the samples of the receive signal are processed in blocks using FFT and 
1 IFFT algorithms. 

15. The method of daim 9 wherein the frequencies, amplitudes and phases of said sinsoidal functions tff time 
are estimated by performing an FFT on the sets of phases of a plurality of pilot tones. 

16. Apparatus for estimating, and compensating for phase jitter imposed on a received multicanier signal by 
imperfect components in a transmission medium, said apparatus comprising: means for multiplying said 
multicanier signal received during a time interval by a locally generated correcting signal to produce a cor- 
rected received multicanier signal, wherein said locally generated correcting signal is an estimate of the 
phase shift required to correct for the effects of said phase jitter 

means for demodulating and then decoding said corrected multicanier signal, estimating the trans- 
mitted data, and then remocJulating said data to form said reference signal; 

means for calculating a phase-error signal which is the difference between the phase of said cor- 
rected multicarier signal and that of said reference signal; 

means for filtering said phase-error signal to generate an estimate of the changes required to the 
parameters of said locally generated correcting signal; 

means for updating said parameters by said changes to generate new parameters of said correcting 
signal; and 

means for repeating the operation of the multiplying, calculating, filtering, and updating means in 
each subsequent time interval. 

1 7. The apparatus of daim 1 6 further including: 

means for sampling said received multicanier signal in synchronism with a locally generated dock 

signal; 

buffer means for assembling the samples of the sampling means into a block; 
means for performing the updating of said parameters after the reception of each block; and 
means for multiplying the samples in a subsequent block of said received multicanier signal by said 
correcting signal, defined by said new parameters. 

18. The apparatus of 17 further including: 

means for calculating said correcting signal by demodulating and then decoding said corrected mul- 
ticanier signal to estimate what data was modulated onto each of the many earners; 

means for filtering the results from said demodulating and decoding means in order to estimate the 
amplitude and corrected phase responses of said transmission medium at each of the carrier frequencies; 

means for remodulating said carriers with said estimated data, using said amplitude and phase res- 
ponses, to generate, as a reference, an ideal received multicarrier signal which has no phase jitter, and 

means for comparing said corrected received multicarrier signal with said remodulated, reference 
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signal to generate said phase-error signal. 

19. The apparatus of claim 18 adapted to the situation in which said correcting signal comprises components 
at several frequencies, called jitter frequencies, at which it has been determined that there is significant 
phase jitter, said apparatus further including: 

means for correlating said phase-error signal with complex sinusoidal signals at said jitter fre- 
quencies; 

means for generating samples of cosine and sine waves of defined amplitudes and phases at each 
of said jitter frequencies; 

means for using the results of said correlation operations to define the changes that should be made 
to said amplitudes and phases; and 

means for summing said samples of cosine and sine waves to generate said correcting signal. 

20. The apparatus of claim 19 wherein said cosine and sine wave generators are implemented by look-up 
tables and multipliers, wherein said look-up tables generate components of said correcting signal with the 
required frequency and phase, and said multipliers adjust the amplitudes of said components. 

21. Apparatus for estimating, and compensating for frequency offset imposed on a received multicarrier signal 
by imperfect components in a transmission medium, said apparatus comprising: 

means for multiplying the received multicarrier signal, received during a time interval, by a locally 
generated correcting signal to produce a corrected received multicarrier signal, wherein said locally gen- 
erated correcting signal is an estimate of the phase shift required to correct for the effects of said frequency 
offset; ' 

means for demodulating and then decoding said corrected received multicarrier signal in order to 
estimate what data was modulated onto each of the many carriers; 

means for filtering the results from said demodulating and decoding means in order to estimate the 
corrected phase response of said transmission medium at each of the carrier frequencies; 

means for fitting and correcting phase response of said transmission medium to a linear function of 
frequency; 

means for interpreting the zero-frequency intercept of said linear function as being caused by a dif- 
ference between the phase of said correcting signal and the phase that is common to all carriers of said 
received multicarrier signal, low-pass filtering said zero-frequency intercept, and using the result of said 
filtering to update the zero-frequency component of said correcting signal. 

22. Apparatus for correcting, in a modem that receives, via a transmission medium, data modulated onto many 
separate carriers, for a difference between the frequencies of a locally generated sampling clock and the 
clock in a remote, transmitting modem, said apparatus comprising: 

means for sampling said received multicarrier signal in synchronism with a locally generated clock 

signal; 

means for interpolating between the resultant samples of said sampling means according to an inter- 
polation parameter, in order to convert samples taken at the frequency of said locally generated clock into 
interpolated samples at the frequency of the clock in said remote transmitting modem; 

means for demodulating and then decoding said interpolated samples in order to estimate what data 
was modulated onto each of the many carriers; 

means for filtering the results of said demodulating and decoding means in order to estimate the 
phase response of said transmission medium at each of the earner frequencies; 

means for fitting said phase response of said transmission medium to a linear function of frequency; 

means for interpreting the slope of said linear function of frequency as being caused by a difference 
between the sampling frequencies of transmitter and receiver; and 

means for low-pass filtering said slope, and using the result as said interpolation parameter. 

23. Apparatus for correcting, in a modem that receives, via a transmission medium, data modulated onto many 
separate carriers, for a difference between the frequencies of a locally generated sampling clock and the 
clock in a remote, transmitting modem, said apparatus comprising: 

means for sampling said received multicarrier signal in synchronism with a locally generated dock 

signal; 

means for controlling the phase of said locally generated dock signal in order to synchronize it with 
the clock in said remote transmitting modem; 
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means for demodulating and then decoding said samples in order to estimate what data was mod- 
ulated onto each of the many carriers; 

means for filtering the results of said demodulating and decoding means in order to estimate the 
phase response of said transmission medium at each of the carrier frequencies; 
5 means for fitting said phase response of said transmission medium to a linear function of frequency; 

means for interpreting the slope of said linear function of frequency as being caused by a difference 
between sampling frequencies of transmitter and receiver; and 

means for low-pass filtering said slope, and using the result to control the phase of said locally gen- 
erated clock. 

10 

24. Apparatus for calculating, in a modem that receives a multicarrier data signal which may be impaired by 
phase jitter, frequency offset and a difference between the frequencies of a locally generated sampling 
clock and the clock in a remote transmitting modem, the initial parameters of a locally generated correcting 
signal which corrects for said impairments, said apparatus comprising: 

15 means for calculating the initial values of the parameters of said correcting signal during the recep- 

tion of a training signal comprising a plurality of pilot tones; 

means for separating said pilot tones by bandpass filters, and calculating their phases at successive 
sampling times; 

means for analyzing said phase to identify linear and sinusoidal functions of time; 
2Q means for estimating the frequencies, amplitudes and phases of said sinusoidal functions of time; 

means for using the estimates produced by said estimating means to initialize the parameters of 
sinewave generators; 

means for comparing the slopes of said linear functions calculated by said analyzing means, and 
for generating an estimate of both said frequency offset and said clock frequency difference; 
25 means for using said estimate of said frequency offset to initialize the zero-frequency component 

of said correcting signal; 

means for using said estimate of the clock frequency difference to initialize the state of a low-pass 
filter and 

means for using the output signal of said low-pass filter to control an interpolator which interpolates 
30 between samples taken in synchronism with a locally generated clock signal to generate samples at the 

frequency of the clock in a remote transmitting modem. 

25. The apparatus of claim 24 wherein the output signal of said low-pass filter is used to control a phase-locked 
loop which adjusts the frequency and phase of said locally generated clock to match that of the clock in a 

35 remote transmitting modem. 

26. The apparatus of claim 16 further including means for processing samples of said received multicarrier sig- 
nal in blocks using Fast Fourier Transform (FFT) and Inverse FFT (IFFT) algorithms. 

40 27. The apparatus of claim 21 further including means for processing samples of said received multicarrier sig- 
nal in blocks using FFT and IFFT algorithms. 

28. The apparatus of claim 22 further including means for processing samples of said received multicarrier sig- 
nal in blocks using FFT and IFFT algorithms. 

45 

29. The apparatus of claim 24 further including means for processing samples of said received multicarrier sig- 
nal in blocks using FFT and IFFT algorithms. 

30. The apparatus of claim 15 further including means for estimating the frequencies of said components of 
so phase jitter by performing an FFT on the sets of phases of said pilot tones. 
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